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On Stability of the Linearized Spacecraft Attitude Control System 


Bin Zhou* 


Abstract 

This note is concerned with the stability and stabilization of the linearized spacecraft attitude control 
system. Necessary and sufficient conditions are respectively provided to guarantee that the considered 
systems are polynomially stable and stable in the Lyapunov sense. These two classes of conditions 
guarantee that the linearized spacecraft attitude control system can be respectively stabilized semi- 
globally and globally by saturated linear state feedback. 

Keywords: Spacecraft Attitude control; Stability and stabilization; Lyapunov stable; Saturated 
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1 Introduction 


Attitude stabilization of spacecraft is a very important problem since it helps the spacecraft to maintain a 
certain prescribed attitude in the space during their useful life n, 0 , a, 0 , 0 , 0 . Attitude stabilization 
controllers design can be based on either the nonlinear dynamics or the linearized dynamics. If the controllers 
are designed based on the nonlinear model, the design procedure can be quite complex while the claimed 
(global) stability can be ensured. If the controllers are designed based on the linear model, the design 
procedure can be much more simpler since a lot of linear design techniques are available, however, the 
claimed (global) stability of the practical closed-loop system may not be maintained. As usual, a trade-off 
thus exists. 

If we are interested in the design of spacecraft attitude stabilizing controllers based on the linearized model, 
it would be very helpful if we can know some specific properties of the linearized model. This is particularly 
the case if we are considering the constraints existing in the actuators since the existence of a controller 
is highly dependent on the stability properties of the open-loop system (see, for example, 13: E: ED: and 
el Hence, in this note, we are interested in the stability and stabilization of the linearized spacecraft 
attitude control system. Particularly, we will provide conditions on the spacecraft parameters under which 
the linearized model is respectively polynomially stable and is stable in the Lyapunov sense. We notice that in 
these two cases the linearized spacecraft attitude control system can be respectively stabilized semi-globally 
and globally by saturated linear feedback (see El and EL 


2 Main Results 


We consider the following linear system 

X = 4\' + Bu , (1) 
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This linear system is the linearized model of the spacecraft attitude control system. For detailed derivation 
of this equation, see a, 0,0, m, and the references therein. In the above equation, the symbol uq = \fl^ 
denotes the orbital rate, in which p, = 3.986 x 10 14 m 3 /s 2 is the gravity constant and r is the semimajor axis 
of the orbit, J x , J y and J z are the inertia of the spacecraft, x contains the attitude angles and velocities, u 
is the forces applied on the spacecraft, and 
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( 3 ) 


In this note, we are interested in the stability properties of system ©■ 


Remark 1 We notice that &i,i — 1,2,3 are dependent. In fact, by letting -f- = fii , -j- = @ 2 , then 
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which implies that ai,i = 1,2,3 can be characterized by the two scalars f3i,i = 1,2. 


Consider the following nonsingular matrices 
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= Wq-Bq, 


where (Ai, A 2 , 2 I 3 ) and (Bx,B 2 , B 3 ) are independent of u >0 and are given by 


( 5 ) 
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Our first result provides necessary and sufficient conditions to guarantee that A is polynomially stable, 
namely, Re{s} < 0,Vs S A (A), the eigenvalue set of A. 


2 






























Proposition 1 All the eigenvalues of A have non-positive real parts if and only if 


( 0 < cr 2 , 

I 0 < (J\a 3 = (pi, 

I 0 < 3(Ji + (J3tTi + 1 = (j>2, 

[ 0 < (3cti + 1 + a 3 ai) 2 — 16cria 3 = (p 2 — 16^i. 

Moreover, if all the above conditions are satisfied, the eigenvalues Si,i G I [1,6] of A are given by 


Si,2 = ±V3o2UJ 0 i, 


S3,4 
S5,6 


<^2 + \/^ 1 - 160 i 

fa -1601 


W 0 1, 

w 0 i. 


(9) 


Proof. Clearly, any s € A (H 2 ) has non-positive real parts if and only if o 2 > 0. On the other hand, if 
s £ A (A 13 ), then it follows from 


det (s /4 — Ai 3 ) = det 


s / 2 -Ax 
—A 3 sh 


det (s 2 / 2 - AxA 3 ) , 


( 10 ) 


that — s € A (A 13 ). Hence Re{s} < 0 if and only if Re{s} = 0. This is equivalent to that all the eigenvalues 
of AiA 3 are non-positive real numbers. Let 

/ (A) = det (A / 2 — A 3 A 3 ) = A 2 + <p 2 X + 4<^>i. (11) 

Hence we must have / (0) > 0, — 4f < 0, A = <f>\ — Idxpi > 0, which are equivalent to ©. Finally, if © are 
satisfied, the eigenvalues of A can be easily computed by noting (fTIIl) . The proof is finished. ■ 

Under the condition that all the eigenvalues of A have non-positive real parts, the linear system (© can be 
semi-globally stabilized by linear feedback in the presence of actuator saturation (see, for example, [2] and 
[TO] ). Particularly, explicit semi-globally stabilizing controllers can be designed by the parametric Lyapunov 
equation approach in HU- 

We next provide necessary and sufficient conditions to guarantee that A is Lyapunov stable (or neutrally 
stable), namely, all the eigenvalues of A have non-positive real parts and those eigenvalues on the imaginary 
axis are simple (or equivalently, any eigenvalue on the imaginary axis has the same algebraic and geometric 
multiplicities). 


Theorem 1 Let A be given by (OH). Then A is Lyapunov stable if and only if (a 1 , er 2 , 03 ) satisfies 

( 0 < (J 2 , 

I 0 < 0i, 

] 0 < <j)2, 

[ 0 < 01 — 16<^i. 


( 12 ) 


Proof. Clearly, A is Lyapunov stable if and only if A 2 and A 33 are all Lyapunov stable. Moreover, the 
matrix A 2 is Lyapunov stable if and only if cr 2 > 0 since A 2 has repeated eigenvalues 0 and is nonzero when 

<J 2 = 0 . 

Obviously, if <f>i > 0 ,i = 1,2,3, then A\ 2 has four different eigenvalues on the imaginary axis and is thus 
Lyapunov stable. Hence, we only need to show the converse, namely, if A 33 is Lyapunov stable, then 
(pi > 0, i = 1,2,3. 

We first show that the Lyapunov stability of A \3 implies (pi > 0. We show this by contradiction. Let 
(pi = 0. If oi = <j 3 = 0, then we must have a 2 = 0 (see equation ©), which is not allowed. Hence we have 
either (ai = 0,a 3 ^ 0) or (ct 3 = 0, cy ^ 0). If (cy = 0, er 3 ^ 0), then A i3 has repeated eigenvalues 0 with 
geometric multiplicity 1, and hence is not Lyapunov stable. If (03 = 0, oy ^ 0), we consider two cases. Case 
1: a 1 ^ —A. Then the algebraic and geometric multiplicities of 0 are respectively 2 and 1, which implies 
that A 33 is not Lyapunov stable. Case 2: a-\ = — A. Then the algebraic and geometric multiplicities of 0 are 
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respectively 4 and 1, which, again, implies that A 13 is not Lyapunov stable. In conclusion, we must have 

<pi > 0 . 


We next show that the Lyapunov stability of A 13 implies (f >2 > 0. Otherwise, if </>2 = 0, then (j > 2 — 16</>i > 0 
is equivalent to </>i < 0, which contradicts with </>i > 0. Hence we must have f >2 > 0. 


We finally show that the Lyapunov stability of A 13 implies (f 2 — 16</>i > 0. Otherwise, if <f\ — 16^>i = 0, then 
A 13 has repeated eigenvalues A 34 = ± \Jfa/^i. Notice that, for any A = Xi ^ 0, * = 3,4, we have 
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where we have noticed that the matrix 

X 2 12 — A 3 A 1 = 


(I - \ a z) ai ~~\ o-i-l 
4 ( 0-3 - 1) fJl \ + (i(T 3 - |) fJl J 7 


(13) 

(14) 


cannot be zero. Hence the geometric multiplicities of A i,i = 3,4, is 1 and A 13 cannot be Lyapunov stable. 
The proof is finished. ■ 


Consider the following Lyapunov equation 

A T P + PA = -D T D , (15) 

where D is any matrix with proper dimensions. It is well known that the above Lyapunov equation has a 
solution P > 0 for some D if and only if A is neutrally stable (Lyapunov stable). Then we have the following 
corollary. 


Corollary 1 There exists a P > 0 such that the Lyapunov equation GW is satisfied for some D if and only 
if D = 0 and \12 1) are satisfied. 

Proof. It follows from Theorem [T| that A is neutrally stable if and only if (fT2l) are satisfied and, in this case, 
all the eigenvalues of A are on the imaginary axis. We next show that D must be zero. Let yt,i £ I [1, 6 ] be 
the eigenvectors of A corresponding to the eigenvalue Sj. Assume that D is not zero. Then there exists a yj 
such that Dyj ^ 0 (otherwise, we must have D = 0). Consequently, it follows from (fl5l) that 

- \\D yj f = -yfD T Dy 3 

= yf(A T P + PA) yj 
= 2 Re (sj)yfPyj 

= 0, (16) 
which contradicts with Dy 3 7 ^ 0. The proof is finished. ■ 

We can actually give explicit solutions to the Lyapunov equation m, as shown in the following theorem. 


Theorem 2 Assume that <j\(T 2 <X 3 7 ^ 0. Then all the solutions to the Lyapunov equation m with D = 0 are 
given by 

" P 2 


P = H 1 P 0 H , P 0 = 
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where cti,a !3 and ai 3 are any scalars and such that 
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Particularly, if we choose ai 3 = 0 and oq = ^pOt. 3 , then we have 
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Pi = diag < a 3 a 3 , —a 3 > , P 3 = diag 4 —CTia 3 , a 3 ^ . 
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Proof. We only need to show that all the solutions Po to the equation Aq Po + PqAq = 0 are in the form 
of (ED. Since Ai,i = 1 , 2 , 3 , are all nonsingular, by noting the structure of A 0 , we know that Po takes the 
following form 

" P 2 0 0 

Pn = 0 Pi P13 , (21) 

0 P13 P3 


which gives the following five equations 


0 — AjP 2 + P 2 A 2 , 

0 = AjP ^ 3 + Pi 3 A 3 , 
0 = Aj P 13 + Pi 3 A\, 
0 = Aj P 3 + P\A\, 

0 = Aj Pi + P 3 A 3 . 


( 22 ) 


We first consider the first equation. Let 


P 2 = 


P2\ P22 

P22 P23 


( 23 ) 


where the parameters are to be specified. Then 


A 2 P2 + P2A2 = 


~6a 2 P22 P21 — 3&2P23 

P21 — 3(J2P23 2p 2 2 


( 24 ) 


Hence we must have p 2 2 = 0 and P21 = 3 ct 2P23■ This is just the desired result if we set P2 3 = 02-with 02 
being any constants. 


We assume that P13 takes the form 


Pi .3 = 


Pl31 P 132 
Pl33 Pl34 


where the parameters are to be specified. Then we can compute 


( 25 ) 


A-J -^13 + P\3A 3 
Aj Pi 3 + P? 3 Ai 


— 2<J 3 Pl 3 2 Pl31 + (o"3 — 1) Pl32 ~ CT3P134 

Pl31 + (&3 — 1) Pl32 — 03P134 2pi33 + 2 (cr 3 — l)pi34 

8(J 1P133 P131 + (1 - CTl)Pl33 - 4(JiPi34 

P131 + (1 - CTi)pi33 - 4aiPi34 2pi32 + 2 (1 - cri)pi34 


( 26 ) 

( 27 ) 


Hence we must have P132 = 0,pi33 = 0 and P131 = 03P134. Consequently, we get 


Aj P\3 + -^13 A 3 
AjPis + P? 3 Ai 


0 0 
0 2 (<r 3 — 1 ) 


P134, 


0 <73 — 4 ai 

OS- 4 ( 7 ! 2 ( 1 -CTl) 


P134- 


( 28 ) 

( 29 ) 


We consider two cases. Case 1 : cti 1. In this case we must have P134 = 0 . Case 2 : cti = 1. Then we have 
(1 — CT3) P134 = 0 and (CT3 — 4)pi34 = 0 . Since ct 3 can not equal to 1 and 4 simultaneously, we also must have 
P134 = 0 . As a result, we have P13 = 0 . 


5 




















We next let Pi and P 3 take the form 


Pll 

P 12 
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where the parameters are to be specified. Then 


AjP, + P[ A-[ 


AJP1 + P3A3 


— C3P32 — 4 (TiPi 2 —CT3P33 + Pn + P12 — CT1P12 

P32 + 0- 3 p 3 2 - 4pi 3 (7i p 3 2 - p 3 3 + &3P33 + Pl2 + Pl3 - P 13 CT 1 J ’ 

— &3P32 — 4<Tipi2 P31 — P32 + CT3P32 — 4pi3(Ti 

—&3P33 + Pll + Pl 2 — CT 1 P 12 P32 — P33 + &3P33 + Pl 2 + Pl3 — P 13 CT 1 


(30) 


(31) 

(32) 


Hence we must have 

{ P 12 = - 03 P 32 / 4 / 0 - 1 , 

Pii = -(- 03 P 33 + P 12 - tripi 2 ), (33) 

P31 = ( P32 + &3P32 - 4pi 3 CTi). 

Consequently, the above two equations are satisfied if and only if 



(1 - G\ ) P 13 + P 32 + (ct 3 - 1) P 33 = 0, 

(34) 

which is 

(1 - (Ji)pi3 + (4(7 1 - 0 - 3 ) P 32 + ( 0-3 - l)p 3 3 = 0. 

(35) 

by setting p 32 4p 32 cri. 

This is just (fl^l) since 



03 - 1 = (ci - 1 ) • 

jZ 

(36) 


and we have denoted P 13 = oq,P 33 = 03 and P 32 = 0 : 12 - The proof is finished by noting that Pi and P3 are 
now exactly in the form of (TTHll . ■ 


Of course, the parameters a\, 0.2 , 013 and <213 can be well chosen such that Pi > 0, * = 1, 2,3, if the conditions 
in m are satisfied. 

The positive definite solutions to the Lyapunov equation m can be utilized to stabilize system m globally 
in the presence if actuator saturation even when B ( t ) is a time-varying periodic matrix. Interested readers 
may refer to [ 12 ] for details. 


3 Conclusion 


This note has studied the stability properties of the linearized spacecraft attitude control system. Necessary 
and sufficient conditions are respectively provided to guarantee that the considered systems are Lyapunov 
stable (neutrally stable) and polynomially stable. 
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